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For a long time, it has been known that the power spectrum of Barkhausen noise had a power-law 
decay at high frequencies. Up to now, the theoretical predictions for this decay have been incorrect, 
or have only applied to a small set of models. In this paper, we describe a careful derivation of 
the power spectrum exponent in avalanche models, and in particular, in variations of the zero- 
temperature random-field Ising model. We find that the naive exponent, (3 — T)/cyvz, which has 
^ ' been derived in several other papers, is in general incorrect for small r, when large avalanches are 

, common, (r is the exponent describing the distribution of avalanche sizes, and avz is the exponent 

describing the relationship between avalanche size and avalanche duration.) We find that for a large 
' class of avalanche models, including several models of Barkhausen noise, the correct exponent for 

r < 2 is 1/avz. We explicitly derive the mean-field exponent of 2. In the process, we calculate the 
(-H ' ■ average avalanche shape for avalanches of fixed duration and scaling forms for a number of physical 

' properties. 
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I. INTRODUCTION 



Many physical systems, from superconductor Jil to sand-pilesi to martensitic shape-memory alloy produce noise 
with power law characteristics. The noise in many of these systems can be modeled in terms of avalanches. One 
example of this is Barkhausen noise in ferromagnetic materials. Many ferromagnetic materials magnetize not smoothly, 
^ . but in jumps of all sizes. The resulting "noise" is characterized by power laws. For example, the distributions of 
Q ■ avalanche sizes, durations and energies are all seen-to-bcLiJower laws. Recently, it has been procpaed that this noise is 
O \ a result of either a disorder induced critical pointclErQiZl'ElEI, or self organized criticality (SOC)IIj'U. Several variations 
of the zero temperature random field Ising model have been proposed to model this critical behavior. 
' One of the main power laws which must be explained if these models are to be successful is the power law behavior 
^ i of the power spectrum. Actually, the power spectrum exhibits two different power laws, one for low frequencies, and 
C another for high frequencies. The high frequency power law, 7^av(w), reflects the dynamics within avalanches, and the 

low frequency power law, 7'corr('^), reflects the correlations between avalanches. Simulations of various random field 
Ising models have been fairly successful in modeling the buzh-frequency scaling of the power spectrum, but theoretical 
^ \ predictions of this scaling have been absent or wrong .ffliJ-EjO (See the discussion in section ^ for a description 
Qv^ ' of several previous calculations.) In this paper, we will derive an exponent relation for the high frequency power 
' spectrum which applies to several variations of the random-field Ising model. Large portions of the derivation should 
also apply to any critical avalanche model. 



II. THE MODELS 



Several variations of the zero temperature random field Ising model have been proposed to explain the power laws 
Q , in Barkhausen noise. They are differentiated on the basis of the presence of long range fopes, and the details of the 
;V ' dynamics. In a separate paper, we examine in detail the differences between these modelslij. A general Hamiltonian 
. !^ , for the models is 



_ J 3cos(%) - 1 

+ / . 2^ -^dipole 3 SiSj, [L) 

where Si — ±1 is an Ising spin, Jnn is the strength of the ferromagnetic nearest neighbor interactions, H^an external 
magnetic field, hi is a random local field, Jinf is the strength of the infinite range demagnetizing field,ll3 and Jdipoie 
is the strength of the dipole-dipole interactions. The critical exponents of the power laws are independent of the 
particular choice of random field distributions p{hi) for a large variety of distributions. Most commonly, a Gaussian 
distribution of random fields is used, with a standard deviation R. (When we refer to the strength of the disorder, 
we are referring to the width, i?, of the random field distribution.) 
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Short Range 


(3D) Front Propagation (3D) Mean Field 


r 


1.6 


1.28 1.5 


avz 


0.58 


0.58 0.5 




Short Range 


(4D) Front Propagation (4D) 


T 


1.53 


1.42 


avz 


0.52 


0.56 



TABLE I: Important exponents for the three universality classes, r is the exponent for the avalanche size distribution D[S) ~ 
avz relates the avalanche size S to the avalanche duration T : T ^ c^avz 



Two different dynamics have been considered. The first is a front propagation dynamics in which spins on the edge 
of an existing front flip as soon as it would decrease their energy to do so. Spins with no flipped neighbors cannot flip 
even if it would be energetically favorable. Second is a dynamics which includes domain nucleation. Any spin can flip 
when it becomes energetically favorable to do so. In both cases, spins flip in shells — all spins which can flip at time t 
flip, then all of their newly flippable neighbors flip at time t + 1. 

Depending on which terms are included in the Hamiltonian, the behavior appears to fall into three different uni- 
versality classes. When domain nucleation is allowed, and only nearest neighbor interactions are included, there is a 
second order critical point at a particular disorder, i?^, and external fleld He- For disorders below R^, a flnite fraction 
of the spins in the Si*t|^W|-teven in the thermodynamic limit) flip in a single jump. This critical point seems to have a 
wide critical region jQEHj'DEtlI so this critical behavior might be seen in experiments even without tuning the disorder. 
For |dwarders below i?c, or when domain nucleation is not allowed, the addition of an inflnite-range demagnetizing 
fleldllSj can self-organize the system to a different critical behavior. (Self-organization means that the system nat- 
urally sits at a critical point, without having to tune any parameters.) The infinite-ranged interaction is sopietimes 
introduced to mimic the effects of the boundaries of materials-with dipolar, or other long-ranged interactionsM; these 
interactions also self-organize the model to the critical pointll3. Our use of the term "inflnite-range models" perhaps 
obscures the clear physical origin of this universality class of models. The criticaLbahavior in front-^jropagation was 
originally described in a non-self-organized depinning model by Ji and Robbins.EZl'EEl Zapperi et aZO claim that the 
addition of dipole-dipole interactions to the inflnite-range model lowers the upper critical dimension to three and 
produces mean-field exponents in three dimensions. Since large mean-field simulations are much easier than large 
simulations with dipole-dipole interactions, we will give results from|^ean-field simulations in this paper. Dipolar 
interactions without an infinite-range term were explored by Magnit2l in two dimensions, who found labyrinthine 
patterns and hysteresis loops similar to those seen in garnet films. The relevant exponents for the three universality 
classes can be found in table ||. Except in section IV, all of the results in this paper are from three-dimensional 
simulations, disorder R = 1.8, nearest neighbor interaction Jnn — 1, and infinite range interaction Ji„f — 0.25; this 
model exhibits self-organized front-propagation exponents. 



III. DERIVING THE EXPONENTS 



In the process of deriving the form of the critical exponent for the power spectrum or energy spectrum (power 
output during the duration of the experiment), wc will derive the scaling forms of several other quantities which are 
themselves of interest. Warned by the failure of the naive scaling exponent for the power spectrum, we will present 
numerical scaling plots for each of these intermediate quantities. 



A. The avalanche shape 



Near criticality, the avalanches in the random-field Ising model have a very ragged shape. There are avalanches 
of all sizes precisely because each avalanche is always finely balanced between continuing and dying out. Most large 
avalanches rGome close to dying many times. The time series for a typical large avalanche can be seen in figure 
Petta et aZB3 have emphasized that the non-trivial structure within avalanches like these are crucial to understanding 
the power spectra, and a good test of the models. 

Despite this rough shape, criticality implies that the average avalanche shape will scale in a universal way. Consider 
the average shape of avalanches of duration T . If we rescale the time axis, i, by a factor of T, and divide the vertical 
axis, which measures the number of spins fiipped (proportionate, in simple cases, to the voltage V which would be 
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FIG. 1; The shape of a typical large avalanche. Notice that the avalanche nearly stopped several times, and the voltage (the 
number of spins flipped in each time step) fluctuated drastically. 

measured in a pickup coil), by the average voltage, we should get a generic shape which is independent of T. The 
average height is the average area, S{T) ~ j^i/o^i^^ this define the exponent avz)f3 divided by the duration, T. 
Therefore, the scaling form should be 

v^(r,t) = Ti/— i/,i,,p,(^/r). (2) 



The scaling of the average avalanche shape according to equation can be seen in figure 




B. Distribution of voltages V in avalanches of size 5* 

The power spectrum is sensitive not only to the shapes of the avalanches, but to the fluctuations in the avalanche 
shapes. An interesting measure of these fluctuations is the probability P{V\S) that a voltage V will occur at some 
point in an avalanche of size S. If this probability scales universally, then we know it must have the form 

p{v\s) = v-^Uo^„{vs-y). (3) 

But what are the exponents x and y? We can determine the value of x by integrating over all voltages. Since P{y\S) 
is a probability distribution, it must integrate to 1: 

v-^Uo^,,{vs-y)dv = As-y^'^-^^ = i. (4) 

From this, we know that x — 1. (The alternative, y = 0, can be discarded because along with equation |^ it would 
imply that P{V\S) is independent of the avalanche size S.) 

We also know that the average voltage in the avalanche must be equal to the avalanche size S divided by the 
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FIG. 2: The average avalanche shape for three different durations. Spasojevic et a/El measured the average avalanche shape 
experimentally and found a somewhat different shape. This kind of measurement provides a much sharper test for the theory 
than the tradition of comparing critical exponents. Presumably, the average avalanche shape for large sizes and times is a 
universal scaling function /shape: if the experiment differs in this regard from our model, then our model is expected to have 
different critical exponents. All features at long length and time scales should be universal. 



avalanche duration T ^ S"''^^ , so 



/"OO 

{V) = / VP{V\S)dV 
Jq 

fpo^cAvs-y)dv 



- Sy - 3^-^"^ (5) 
From this, we know that y = I — avz, and the probability of a voltage V occurring in an avalanche of size S is 

p{y\S) = V^-Vpower(^^"''^-'). (6) 



As can be seen in figure ^ this scaling form works very well. 



C. The scaling of avalanche energy E with avalanche size 5* 

The voltage distribution in equation |^ allows us to calculate the dependence of avalanche energy on avalanche size. 
The avalanche energy is simply the average squared voltage, {V"^), times the average avalanche duration, 5°^"^. Using 
equation ^ this is 



/"OO 

E{S) = S""'' / V^P{V\S)dV 
Jo 



/•OO 

S'"'' / Vfpo^,r{VS'''''-^)dV 

Jo 

c2 — erf 2 



(7) 
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FIG. 3: A collapse of the voltage distribution for three avalanche sizes according to equation ^ 
Note that this is the same result we would find if we assumed that the time dependence had a square profile. 

D. The scaling of the time-time correlation function with S 

With this information, we can calculate the scaling of the time-time correlation function within avalanches, which 
is simply related to the high-frequency part of the power-spectrum, Pav(^)- The time-time correlation function is 
defined as 

G(B) = J V{t)V{t + 9)dt. (8) 

If we assume that the magnetic field is increased adiabatically and the avalanches are well separated in time, we 
can calculate the time-time correlation function separatelA2-|for each avalanche and then add the individual functions 
together to get the overall time-time correlation function.c3 

This allows us to break up the time-time correlation function into the contributions from avalanches of different 
sizes, S. Let G{0\S) be the average time-time correlation function of an avalanche, given that the avalanche is of size 
S. (In contrast, the notation G{9,S) would denote the contribution of all avalanches of size S to G{9). This would 
be weighted by the probability that an avalanche is of size S*, S~'^ .) If we consider equation || at 6* = 0, we see that 
the 9 — Q component of the correlation function is J V(ty'dt, which is proportional to the avalanche energy. Using 
this fact along with the scaling of the energy from equation ^ we find that the time-time correlation function should 
scale as 

G{9\S) = S^-^'''hr..corAOS-^n- (9) 

As shown in figure ^, this scaling works very well over a wide range of avalanche sizes. 

E. The scaling of the power spectrum with S 

The power or energy spectrum describes the amount of energy released in Barkhausen noise at each frequency. It 
can be calculated as the cosine transform of the time-time correlation function. Transforming equation ^ we find that 
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FIG. 4: A collapse of avalanche time-time correlation functions according to equation ||. The data is binned logarithmically to 
get good statistics. 



the scaling of the power spectrum with avalanche size has the form 

/>oo 

E{uj\S) = / cos{uj9)G{9\S)de 



COs{cO0)S^~^'''ftuaccoA0S-''''')d0 

_ q2 r ( l/(7UZ Q\ 

— ^ ^timccorrV'-^ 



(10) 



(We'll call this the power spectrmii, but use the name E[ijj) to remind us to divide by the duration of the experiment 
to change energy to power.) A collapse of the power spectra for different avalanche sizes according to this form is 
shown in figure pi 



F. The power spectrum scaling depends on the large avalanche cutoff 

The power spectrum for all avalanches (the quantity usually measured in experiments) is just the integral 
/ D{S)E{Ljj\S)dS , where D{S) ~ S~'^ is the probability per unit volume of having an avalanche of size S. If we 
substitute equation and D(S) ~ S^'^ into the integral, we generate the naive prediction of 

i?wrong(c^) = f DiS)E{u;\S)dS ^ Cc>-(3-r)/a.. [ y^- f,^,,^{y)dy ^ c^'O-)/-- (11) 



as predicted by Lieneweg et aL.Ill Dahmen and Sethna,Bi@ and Spasojevic et azJll. In the next section [II G we will 
argue that the energy contributed by an avalanche at a fixed frequency to is proportional to S: this immediately 
implies using the scaling form |l^ that /cnorgy(2/) ~ A/y for large y. Thus the integrand in equation |l^ dies off as a 
power law y^^'^ for large y. Hence, it gives the correct answer if r > 2, but for t < 2 the naive indefinite integral in pi] 
diverges and must be replaced by a definite integral. The naive result should work for the integrated avalanche size 
distribution for the short-range modeloCI at the critical point R^- In that case (which we won't study in this paper), 
the corresponding exponent t + al35 — 2.03 in three dimensions (close to two, so we'd expect logarithmic corrections), 
and rising to the mean-field value of 9/4 in six dimensions. For the models we study here, r < 2 and the cutoff at the 
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FIG. 5: A collapse of the power spectra for different avalanche sizes according to equation hol 



largest avalanches in the integral over avalanche sizes changes the form of the pov^fer spectrum. In the section [II H 
after next, we will derive our resulting exponent relation for the power spectrum. 



G. The energy at frequency ui scales linearly with the avalanche size S. 

We shall argue in this section that except at very small frequencies the energy at a fixed frequency lu is proportional 
to S. According to equation |lO|, this implies that the power spectrum at fixed S scales as uj~^^'"^'^ . This means that 
dividing E{lj\S) by Slj^^^'^'^^ should collapse the curves and eliminate the w dependence. The result of this collapse 
for the infinite-range model can be seen in figure ^. 

The fact that the energy at frequency lu scales linearly with S follows from two assumptions. First, each spin in 
the avalanche contributes equally to E{uj). Second, the contribution of a given spin is independent of the size of the 
avalanche it is in. (Of course, a few spins at the beginning of the avalanche might differ, so long as the fraction of such 
spins vanishes for large avalanches.) One way in which this could be true would be if the following two hypotheses 
were true. First, each spin contributes to E[lij) only through its correlations with physically nearby spins. Second, the 
local growth of an avalanche does not reflect the overall size of the avalanche. The first hypothesis seems likely in the 
models described in this paper because the times of physically distant spin flips are likely to be randomly distributed 
in time and contribute incoherently to the power spectrum. Only nearby spin flips will be correlated in time and 
contribute coherently. The second hypothesis is also likely to be true because the avalanches are occurring at a critical 
point. Every part of the avalanche is always on the verge of stopping regardless of the size of the avalanche. 

We can check these hypotheses by breaking up the power spectrum into the contributions from pairs of spins at 
different radii r. The time-time correlation function for avalanches of size S can be written in terms of individual 
pairs of spins as 



G{e\s) 



hi 



(12) 



where ti is the time at which spin i flips. From this form, we can use an additional delta function to pull out the 
contribution due to pairs of spins separated by a distance r: 



G{e,r\S) = (Y,Sit,-U-9)Si\r, 



(13) 
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FIG. 6: A collapse of the power spectra for several avalanche sizes S. The collapse is performed by dividing out the linear S 
dependence. The curves are made flat by dividing out the simple power law a;^'''^''^. The bj axis has not been rescaled, so the 
high uj cutoffs collapse together and the small ci; cutoffs do not. 



where ()i implies an averaging over all values of i. Using this definition of the function G(B^ ''|<5'), we can rewrite the 
time-time correlation function as 

G{B\S) = S J G{e,r\S)dr. (14) 

Now, we can use equation |lj to calculate the contribution of spins separated by a distance r to the power spectrum 
E{uj\S) at a frequency uj. Taking the cosine transform of equation we find that 



E{lu\S) = S J J co8{ioe)G{e,r\S)drde 

= S I E{uj,r\S)dr, (15) 



where the function E{uj, r \S) is defined by this equation. Notice that E{uj, r \S) must have a cutoff at the largest r 
present in an avalanche of size S. We can see from equation |l^ that in order for E(uj\S) to be proportional to S, 
the integral must be independent of this S dependent cutoff. This is a more precise statement of the hypothesis that 
only correlations between nearby spins contribute to the power spectrum. It is also necessary that except for extreme 
values of uj and r, E(uj,r \S) must be independent of S. Combined with the first condition, this corresponds to the 
hypothesis that the local growth of the avalanche should not refiect the overall size of the avalanche. 

Figure ^ shows E{uj,r \S) is decaying approximately as 1/r. Figure || shows that it is independent of S and that 
the decay is oscillating about zero. Neighboring radial shells contribute with opposite sign. Hence the integral |l^ 
appears to be conditionally convergent at large distances: each spin contributes to E{uj) only through its correlations 
with nearby spins, and E{uj\S) ~ S. 



H. Integrating the power spectrum over avalanche sizes 



In section IIIF we found that we need to understand the cutoff at large avalanche sizes. In our simulations there is 
a typical largest avalanche size S* which is primarily due to finite size effects: S* ^ L^l"'^ . In experiments, the typical 
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FIG. 7: This graph shows the function \E{Lj,r \S)\ for a range of values of r at 5 = 32500. The function E{uj,r\S) not only 
decays with uj, but also oscillates about zero. To better compare the amplitudes of the curves at different r, we took the 
absolute value and performed a rumiing average over u, averaging out the oscillations. The horizontal lines show how the 
positions of these curves should scale if they went as 1/r. Notice that while for small r the amplitudes drop slower than 1/r, 
for larger r the amplitudes drop off approximately as 1/r. 
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FIG. 8: This graph shows r \S)\ at r = 6 for two sizes of avalanches. No averaging has been performed, so the oscillations 

are visible (the function changes sign at each dip). For larger r, the oscillations are much faster. Notice that the curves are 
nearly identical for the two sizes — E{u!,r \S) has no significant S dependence at small r. 
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largest avalanche size could be determined by the finite experimental duration or from the demagnetizing forces. (In 
some cases, demagnetizing forces also contribute to the cutoff in our simulations.) These effects would cut off the 
probability of getting large avalanches. 

The contribution to E{lu) from an avalanche of size S scales linearly in S, for large avalanches: each spin contributes 
the same amount. The cutoff at large avalanche sizes dominates the scaling of E{uj) precisely because most spins 
are in the largest avalanches. The distribution of avalanche sizes D(S, S*) cannot be a simple power law cut off at 

S* for T < 2, because we know / SD{S,S*)dS = 1; J" S^-'-dS depends on S* and diver ges as 5* — > oo for t < 2. 
Hence the overall amplitude of D must decreases as S* gets bigger. If we choose a scaling form for the avalanche size 
distribution with an overall amplitude which depends upon S* , we can fix this problem: 

D{S, S*)dS = {S*r^UUS/S*) (16) 

fsi.eiy)-y-^ y^O (17) 

J SD{S, S*)dS = J {S*)-^SU,c{S/S*)dS = J yUUy)dy = 1 (18) 

where the last equation provides a condition on the scaling function /size. 

We are now ready to integrate the power spectrum over S*, finding the form E[uj) corresponding to equation |ll| 
but valid for r < 2. We discovered in the last sections that /encrgy(j/) ^ -^/u- We control the integral by adding and 
subtracting A/y: 

E{uj)= j D{S, S*)E{uj\S)dS 

[{S*)-^UUSIS*)][S^f.n..sy{^^''''''S)]dS 

{S/S*fUUS/S*)foncrsy{^'/"'''S)dS 

{S/S*fUUS/S*)A/{u;^/^'''S)dS 



+ I iS/S*ffsUS/S*)[fener,yiio'/^'''S) - A/ (cj' S)]dS (19) 



In the first integral, we set z = S/S*. In the second integral, we set y = uj^/'^'^^S. Also, the second integral now 
converges, so we can substitute fsizeiS/S*) (S/S*)^'^: 

E{lo) = cc^-i/-- I Az^UU^.)dz + (5*)-2^-(3-r)/a.. f y^-r^j^^^^^^^y) _ Aly]dy (20) 



(t < 2) (21) 

The exponent in the second term of equation (3 — t)/ avz, is -Ewrong of equation given by ignoring the cutoff 
at large avalanche sizes. However, for t < 2 the first term will dominate over the second term both for large system 
sizes and for large oj. Only for r > 2 will the second term dominate. For all of the models we are considering, r < 2. 
(See table |.) In fact, the exponent (3 — T)/avz disagrees badly with the results observed in both simulations and 
experiments, and the exponent l/avz agrees very well. In mean field theory, we can actually derive rigorously that the 
power spectrum exponent is 2 (see appendix |A[). This is in perfect agreement with the exponent l/avz, and complete 
disagreement with the exponent (3 — T)/avz, which would be 3. A plot of the power law for the infinite-range model 
can be seen in figure ^. A comparison of the two possible exponents with simulations for each of the three models can 
be seen in table ||. 



IV. HOW UNIVERSAL IS THE EXPONENT 1/auzl 



So far in this paper, we have only discussed these results for the infinite range model in three dimensions. However, 
we have also checked these results in several other situations where both r and avz take on a range of values. (In 
all cases, r < 2.) We have checked the short-range model in 3 and 4 dimensions, the infinite-range model in 3 and 4 
dimensions, and mean field theory. In all cases, the linear scahng with S seems exact, and both E[ijo\S) and E{uj) scale 
as 0;"^/°''^^ to within simulational precision. In all cases, the exponent (3 — t)I avz is completely inconsistent with the 
observed results. The results of mean field theory and the short-range model at the three dimensional critical point 
can be seen in figures |l^ and |ll]. The results for the four dimensional short-range model and the four dimensional 
infinite range model were also completely consistent with a Xjavz scaling exponent. 
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FIG. 9: The power spectrum for the infinite range model. The dashed line is not a fit. It is a power law with an exponent 
Xjavz. 





Ijavz (3 


~ r) javz 


Simulated 


Short-range 


1.72 


2.41 


1.70 


Infinite-range 


1.72 


3.00 


1.70 


Dipole 


2.00 


3.00 


2.00 



TABLE II: This table compares values of the power spectrum exponent measured in simulations with the values predicted by 
the exponents Xjavz and (3 — T)/auz. Notice that for all three models 1/auz is very close to the simulated exponent, and 
(3 — T^javz is in complete disagreement. 



V. DISCUSSION 

There have been several previous predictions of the power spectrum exponent. In her Ph.D. thesis, Karin DahmenE 
did a calculation similar to the one done in this paper, but she ignored the problems with the integral and came up 
with tiie exponent (3 — t)/ avz, which we have shown should hold only if t > 2. This exporvB«|t-was published by ohs 
groupO without the derivation, and was compared to experimental work by Cote and MeiseOEa and Bertotti et alH 
The correct exponent form derived here makes the agreement between the short-range model and these experiments 
substantially better: they quote an exponent of "around 2", our former (wrong) prediction was 2.46 ± 0.17, and the 
correct prediction is 1.75 ± 0.25. 

Presumahlji the earliest derivation of (3 — T^javz (we thank the referee for pointing this out), is b« Lieneweg and 
Gross-NobisE3 in 1972: their exponent 1 -I- 2e — (5 translates precisely to (3 — T)/avz in our notation. t3 They use this 
to analyze their unannealed data, for which they measure r = 2.1 (hence making their derived exponent appropriate); 
their annealed wires show values of r = 1.73, and indeed their power spectrum shows a much steeper falloff, quite 
close to our predicted value \/uvz which thw, measure to be 1.63. 

Spasojevic, Bukvic, Milosevic and Stanleyli^ also came up with the same (wrong) form based on arguments about 
the average pulse shape. They measured an experimental power spectrum exponent of 1.6. For their experimentally 
determined values of r = 1.77 and avz = 0.662, they found that (3 — T)/avz — 1.86 fit their data well, but the value 
of (3 — T)/avz 2.46 quoted by Dahmen et. aLB was much too large. This lead them them to disregard the plain 
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FIG. 10: The mean field model. Notice that the power spectrum for avalanches of size S collapses when divided by S, and 
that the power spectrum for avalanches of size S and the overall power spectrum both have the same power law as the average 
avalanche size as a function of the inverse avalanche duration. Sizes 96, 1536, 24576, and 196608 are shown in the collapse of 
P{uj\S)/ S. The average avalanche size as a function of the inverse avalanche duration is the definition of the exponent 1/auz. 
Power laws of l/ai/z and (3 — T)/auz are also shown for comparison. Notice that (3 — T)/auz is completely incompatible with 
the results. 




FIG. 11: The short-range model in 3 dimensions. The plots and the sizes are the same as in figure y_0| 
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old critical model with domain nucleation as a possible explanation of Barkhausen noise in their system. However, 
note that their experimental value of X/avz — 1.51 is as close to their measured power spectrum exponent as their 
value of (3 — t)/(7vz ~ 1.86. The value of 1/avz — 1.75 ± 0.25 quoted by Dahmen et al. is even closer to the 
experimentally observed power spectrum exponent. (The value of t = 1.6 predicted by Dahmen et al. is also closer to 
the experimentally measured r — 1.77 than any of the other models: front propagation has an exponent r = 1.3, and 
mean-field theory has an exponent r = 1.5.) One should note, however, that the average avalanche shapes measured 
by Spasojevic et al. disagree with those prec^ted by the short-range model (figure 

In 1989, Jensen, Christensen and FogedbyEZI published a different calculation of the power spectrum exponent for 
sandpile models, which has been cited many times as an explanation for 1/w^ scaling of the power spectrum in 
sand-piles, Barkhausen noise, and other systems. They made two major assumptions in their derivation. First, they 
assumed that the avalanche shape could be approximated by a box function: V{t) ~ S/T for all t < T. For our models, 
this assumption turns out to be valid for calculating the average avalanche energy, but is not valid for determining the 
overall scaling of the time-time correlation function. Second, they assumed that one of their scaling functions, which 
was related to the time-time correlation function, had the simple scal ing form G{T) ~ T" exp(— T/Tb). As can be seen 
from the exact mean field time-time correlation function in equation A6, this is not at all a safe assumption. In terms 
of our exponents t and ai^z, their prediction was that the power spectrum would have the form E{uj) ^ uji^^'^)/'^'^^ 
for r > 3 — 2ai'z and E{uj) ^ uj^'^ for r < 3 — 2(ji'z. This is the same as our result in the case where avz — 1/2. 
For the particular model they were studying, their assumptions seem to be correct, but in the more general case, the 
exponent (3 — T)/avz should become invalid for r < 2, rather than r < 3 — avz. Also, for t < 2, exponents other 
than 2 are possible, depending on the particular dynamics. For the class of models described in this paper, the correct 
exponent for. r < 2 is l/avz. 

NarayanES derives a power-spectrum exponent for fcont propagation, considering the dynamic height- height correla- 
tion function for the growing interface. His exponentiHa, {2(^ + d— l)/z — 1 involves the static height-height correlation 
exponent C: {h{x,t)h' {x' ,t)) ~ jx — x'p'', along with the dimension d. This^rediction is rather different from the two 
discussed in this paper. Narayan and Fisher's 5 — d = e-expansion resultsB2l show ( = 2/3 and z = 2 — 2e/9 + O(e^) 
yielding a predicted power-spectrum exponent of 8/7 -I- O(e^) in three dimensions, quite different from the value of 
around 1.7 that we measure. Narayan's derivation does not separate the high-frequency power law 'Pav(t^) which 
we study from the low-frequency power law 7'corr('^) reflecting the correlations between avalanches (discussed in the 
introduction): his analysis is probably describing the asymptotic behavior of the latter. 

There are other specific models for which the power spectrum exponent has been calculated. For example, Bak, 
Tang and WiesenfeldH calculated the power spectrum exponent for their sandpile models, and—Diirin, Bertotti, and 
Magncj calculate the form and asymptotic power law for the power-spectrum of the ABBMMEj model. Without 
further investigation, we can't expect our arguments for an exponent of l/avz to hold for these and other models. 
However we do expect that in any avalanche based model, the power spectrum exponent for r > 2 will be (3 — r)/ avz, 
and another exponent will dominate for r < 2. Whenever the arguments in section III G hold, for r < 2 the exponent 
l/avz, relating the avalanche size to avalanche duration, should dominate. 
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APPENDIX A: THE MEAN FIELD POWER SPECTRUM 



In mean field theory, we can calculate the power spectrum exactly. The Hamiltonian in mean field theory is the 
Hamiltonian in equation ^ without nearest neighbor or dipole terms: 

H = -Y2iH + J + h^)s^. (Al) 

i 

When a spin flips with an external field of H, all spins with random fields between —{H+JM) and —{H+J{M+2/N)) 
will flip. Therefore, each spin has a probability of ^p{—H — JM) of flipping, where p{h) is the probability distribution 
of the random fields. On average, 2Jp{—H — JM) spins will be flipped. If 2J p{—H — JM) > 1, then the avalanche 
will tend to grow indefinitely, and there will be an infinite avalanche. If 2J p{—H — JM) < 1, then the avalanche 
will quickly die out, and all avalanches will be small. If 2Jp{—H — JM) — 1, then avalanches will always be finely 
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balanced between continuing and dying, and there will be a critical distribution of avalanches. If the random field 
distribution p{h) has a maximum value of 1/2 J, then there will be a critical distribution of avalanches at the value of 
H = He where p{h) is a maximum. 

Let us calculate the power spectrum for the critical system with p{h) ~ 1/2 J. To begin with, we will calculate the 
probability distribution of time-series rti,n2, . . . ,noo- We know that in shell zero, exactly hq — I spins will flip. At 
the critical point, where each spin on average causes one more spin to flip, shell one will have a Poisson distribution 

with mean one: P{ni) = Shell i will have a Poisson distribution with mean ni_i: P{ni) — Therefore, 

the probability distribution for the entire time series will be 



2 °° e"*"in"' 

P(l, n„n„...,n^)^—J[ (A2) 



Now, from equation A2, we can calculate the average time-time correlation function 

oo 

G{e) = (A3) 

i=0 

where 

{n^Ui+g) ^ ^ nin^+gP{ni,...,noo)- (A4) 

{nl,...,n^}=0 

To simplify equation |A^ , we need to use several properties of the Poisson distribution: J2^=a ^ „f = 1, 

X^^o*^^ — r' ~ ^n=o T^r" ~ X + x^. Using the first rule repeatedly, we can simplify equation A_4 

to 

OO 

{niUi+g) ^ ^ nifli+gPim, Tli+g). 

{ni,...,ni+8}=0 

Then, using the second rule repeatedly, we can further simplify to 

oo 

{niUi+g) = ^ nfP{ni,...,ni). 

{ni ,...,ni}— 

Now, applying the second and third rules repeatedly, we can simplify to a single sum: 

{i — 1)77,1 + ni 



(jliTli+g) = ^ 



eni ! 

rii=0 ^ 



= i + l. (A5) 

Notice that the correlation between two times is proportional only to the first time, and not the separation between 
the times. 



Now, we can find the value of the time-time correlation function G{9) by summing the result of equation A5 



Because the time-time correlation function as defined is proportional to the square of the total time, we must cut off 



the summation at some maximum time T to get a finite result. Summing according to equation A3, we find 

i=0 

_ {T -e){T-e + i) 
~ 2 

2222 ^ ' 

This shape for a cutoff time of T = 1000 is shown in figure |l^. Notice that the exact shape in mean-field theory is 
very similar to the experimentally measured correlation functions in three dimensions. 
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FIG. 12: The exact time-time correlation for mean field theory with a maximum time cutoff of T = 1000. 



To get the exact form of the mean-field power spectrum, we just take the cosine transform of equation A6 

E{oj) = / cos{ujd)G{0)de 
Jo 



Tlj- 



(1 - cos(rw)) - uj-^ sin(rw). 



(A7) 



The dominant term in this equation is Tlu^'^. This is the same exponent as predicted by the general scaling arguments: 
-^j^ = 2. The general scaling arguments also predicted that a term smaller by a factor of T and with the exponent 
(3 — t)I<jvz = 3 would be subtracted off, and it is, but multiplied by a factor of sm^TtS). However, the sin(rti;) turns 
out to simplify things even more. Because the correlation function was actually discrete, we should consider a discrete 
power spectrum, where the frequencies are multiples of = This means that cos(Ta;) = 1 and sin(Tci;) = for 
all u) in the discrete spectrum. Because of this, all terms except the Tuj~'^ term drop out. If we divide by T to get 
the power spectrum, we find 



(A8) 
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For a long time, it has been known that the power spectrum of Barkhausen noise had a power-law 
decay at high frequencies. Up to now, the theoretical predictions for this decay have been incorrect, 
or have only applied to a small set of models. In this paper, we describe a careful derivation of 
the power spectrum exponent in avalanche models, and in particular, in variations of the zero- 
temperature random-field Ising model. We find that the naive exponent, (3 — T)/cyvz, which has 
^ ' been derived in several other papers, is in general incorrect for small r, when large avalanches are 

, common, (r is the exponent describing the distribution of avalanche sizes, and avz is the exponent 

describing the relationship between avalanche size and avalanche duration.) We find that for a large 
' class of avalanche models, including several models of Barkhausen noise, the correct exponent for 

r < 2 is 1/avz. We explicitly derive the mean-field exponent of 2. In the process, we calculate the 
(-H ' ■ average avalanche shape for avalanches of fixed duration and scaling forms for a number of physical 

' properties. 



C/3 



a 

X3 



O 



I. INTRODUCTION 



Many physical systems, from superconductor Jil to sand-pilesi to martensitic shape-memory alloy produce noise 
with power law characteristics. The noise in many of these systems can be modeled in terms of avalanches. One 
example of this is Barkhausen noise in ferromagnetic materials. Many ferromagnetic materials magnetize not smoothly, 
^ . but in jumps of all sizes. The resulting "noise" is characterized by power laws. For example, the distributions of 
Q ■ avalanche sizes, durations and energies are all seen-to-bcLiJower laws. Recently, it has been procpaed that this noise is 
O \ a result of either a disorder induced critical pointclErQiZl'ElEI, or self organized criticality (SOC)IIj'U. Several variations 
of the zero temperature random field Ising model have been proposed to model this critical behavior. 
' One of the main power laws which must be explained if these models are to be successful is the power law behavior 
^ i of the power spectrum. Actually, the power spectrum exhibits two different power laws, one for low frequencies, and 
C another for high frequencies. The high frequency power law, 7^av(w), reflects the dynamics within avalanches, and the 

low frequency power law, 7'corr('^), reflects the correlations between avalanches. Simulations of various random field 
Ising models have been fairly successful in modeling the buzh-frequency scaling of the power spectrum, but theoretical 
^ \ predictions of this scaling have been absent or wrong .ffliJ-EjO (See the discussion in section ^ for a description 
Qv^ ' of several previous calculations.) In this paper, we will derive an exponent relation for the high frequency power 
' spectrum which applies to several variations of the random-field Ising model. Large portions of the derivation should 
also apply to any critical avalanche model. 



II. THE MODELS 



Several variations of the zero temperature random field Ising model have been proposed to explain the power laws 
Q , in Barkhausen noise. They are differentiated on the basis of the presence of long range fopes, and the details of the 
;V ' dynamics. In a separate paper, we examine in detail the differences between these modelslij. A general Hamiltonian 
. !^ , for the models is 



_ J 3cos(%) - 1 

+ / . 2^ -^dipole 3 SiSj, [L) 

where Si — ±1 is an Ising spin, Jnn is the strength of the ferromagnetic nearest neighbor interactions, H^an external 
magnetic field, hi is a random local field, Jinf is the strength of the infinite range demagnetizing field,ll3 and Jdipoie 
is the strength of the dipole-dipole interactions. The critical exponents of the power laws are independent of the 
particular choice of random field distributions p{hi) for a large variety of distributions. Most commonly, a Gaussian 
distribution of random fields is used, with a standard deviation R. (When we refer to the strength of the disorder, 
we are referring to the width, i?, of the random field distribution.) 
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Short Range 


(3D) Front Propagation (3D) Mean Field 


r 


1.6 


1.28 1.5 


avz 


0.58 


0.58 0.5 




Short Range 


(4D) Front Propagation (4D) 


T 


1.53 


1.42 


avz 


0.52 


0.56 



TABLE I: Important exponents for the three universality classes, r is the exponent for the avalanche size distribution D[S) ~ 
avz relates the avalanche size S to the avalanche duration T : T ^ c^avz 



Two different dynamics have been considered. The first is a front propagation dynamics in which spins on the edge 
of an existing front flip as soon as it would decrease their energy to do so. Spins with no flipped neighbors cannot flip 
even if it would be energetically favorable. Second is a dynamics which includes domain nucleation. Any spin can flip 
when it becomes energetically favorable to do so. In both cases, spins flip in shells — all spins which can flip at time t 
flip, then all of their newly flippable neighbors flip at time t + 1. 

Depending on which terms are included in the Hamiltonian, the behavior appears to fall into three different uni- 
versality classes. When domain nucleation is allowed, and only nearest neighbor interactions are included, there is a 
second order critical point at a particular disorder, i?^, and external fleld He- For disorders below R^, a flnite fraction 
of the spins in the Si*t|^W|-teven in the thermodynamic limit) flip in a single jump. This critical point seems to have a 
wide critical region jQEHj'DEtlI so this critical behavior might be seen in experiments even without tuning the disorder. 
For |dwarders below i?c, or when domain nucleation is not allowed, the addition of an inflnite-range demagnetizing 
fleldllSj can self-organize the system to a different critical behavior. (Self-organization means that the system nat- 
urally sits at a critical point, without having to tune any parameters.) The infinite-ranged interaction is sopietimes 
introduced to mimic the effects of the boundaries of materials-with dipolar, or other long-ranged interactionsM; these 
interactions also self-organize the model to the critical pointll3. Our use of the term "inflnite-range models" perhaps 
obscures the clear physical origin of this universality class of models. The criticaLbahavior in front-^jropagation was 
originally described in a non-self-organized depinning model by Ji and Robbins.EZl'EEl Zapperi et aZO claim that the 
addition of dipole-dipole interactions to the inflnite-range model lowers the upper critical dimension to three and 
produces mean-field exponents in three dimensions. Since large mean-field simulations are much easier than large 
simulations with dipole-dipole interactions, we will give results from|^ean-field simulations in this paper. Dipolar 
interactions without an infinite-range term were explored by Magnit2l in two dimensions, who found labyrinthine 
patterns and hysteresis loops similar to those seen in garnet films. The relevant exponents for the three universality 
classes can be found in table ||. Except in section IV, all of the results in this paper are from three-dimensional 
simulations, disorder R = 1.8, nearest neighbor interaction Jnn — 1, and infinite range interaction Ji„f — 0.25; this 
model exhibits self-organized front-propagation exponents. 



III. DERIVING THE EXPONENTS 



In the process of deriving the form of the critical exponent for the power spectrum or energy spectrum (power 
output during the duration of the experiment), wc will derive the scaling forms of several other quantities which are 
themselves of interest. Warned by the failure of the naive scaling exponent for the power spectrum, we will present 
numerical scaling plots for each of these intermediate quantities. 



A. The avalanche shape 



Near criticality, the avalanches in the random-field Ising model have a very ragged shape. There are avalanches 
of all sizes precisely because each avalanche is always finely balanced between continuing and dying out. Most large 
avalanches rGome close to dying many times. The time series for a typical large avalanche can be seen in figure 
Petta et aZB3 have emphasized that the non-trivial structure within avalanches like these are crucial to understanding 
the power spectra, and a good test of the models. 

Despite this rough shape, criticality implies that the average avalanche shape will scale in a universal way. Consider 
the average shape of avalanches of duration T . If we rescale the time axis, i, by a factor of T, and divide the vertical 
axis, which measures the number of spins fiipped (proportionate, in simple cases, to the voltage V which would be 
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FIG. 1; The shape of a typical large avalanche. Notice that the avalanche nearly stopped several times, and the voltage (the 
number of spins flipped in each time step) fluctuated drastically. 

measured in a pickup coil), by the average voltage, we should get a generic shape which is independent of T. The 
average height is the average area, S{T) ~ j^i/o^i^^ this define the exponent avz)f3 divided by the duration, T. 
Therefore, the scaling form should be 

v^(r,t) = Ti/— i/,i,,p,(^/r). (2) 



The scaling of the average avalanche shape according to equation can be seen in figure 




B. Distribution of voltages V in avalanches of size 5* 

The power spectrum is sensitive not only to the shapes of the avalanches, but to the fluctuations in the avalanche 
shapes. An interesting measure of these fluctuations is the probability P{V\S) that a voltage V will occur at some 
point in an avalanche of size S. If this probability scales universally, then we know it must have the form 

p{v\s) = v-^Uo^„{vs-y). (3) 

But what are the exponents x and y? We can determine the value of x by integrating over all voltages. Since P{y\S) 
is a probability distribution, it must integrate to 1: 

v-^Uo^,,{vs-y)dv = As-y^'^-^^ = i. (4) 

From this, we know that x — 1. (The alternative, y = 0, can be discarded because along with equation |^ it would 
imply that P{V\S) is independent of the avalanche size S.) 

We also know that the average voltage in the avalanche must be equal to the avalanche size S divided by the 
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FIG. 2: The average avalanche shape for three different durations. Spasojevic et a/El measured the average avalanche shape 
experimentally and found a somewhat different shape. This kind of measurement provides a much sharper test for the theory 
than the tradition of comparing critical exponents. Presumably, the average avalanche shape for large sizes and times is a 
universal scaling function /shape: if the experiment differs in this regard from our model, then our model is expected to have 
different critical exponents. All features at long length and time scales should be universal. 



avalanche duration T ^ S"''^^ , so 



/"OO 

{V) = / VP{V\S)dV 
Jq 

fpo^cAvs-y)dv 



- Sy - 3^-^"^ (5) 
From this, we know that y = I — avz, and the probability of a voltage V occurring in an avalanche of size S is 

p{y\S) = V^-Vpower(^^"''^-'). (6) 



As can be seen in figure ^ this scaling form works very well. 



C. The scaling of avalanche energy E with avalanche size 5* 

The voltage distribution in equation |^ allows us to calculate the dependence of avalanche energy on avalanche size. 
The avalanche energy is simply the average squared voltage, {V"^), times the average avalanche duration, 5°^"^. Using 
equation ^ this is 



/"OO 

E{S) = S""'' / V^P{V\S)dV 
Jo 



/•OO 

S'"'' / Vfpo^,r{VS'''''-^)dV 

Jo 

c2 — erf 2 



(7) 
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FIG. 3: A collapse of the voltage distribution for three avalanche sizes according to equation ^ 
Note that this is the same result we would find if we assumed that the time dependence had a square profile. 

D. The scaling of the time-time correlation function with S 

With this information, we can calculate the scaling of the time-time correlation function within avalanches, which 
is simply related to the high-frequency part of the power-spectrum, Pav(^)- The time-time correlation function is 
defined as 

G(B) = J V{t)V{t + 9)dt. (8) 

If we assume that the magnetic field is increased adiabatically and the avalanches are well separated in time, we 
can calculate the time-time correlation function separatelA2-|for each avalanche and then add the individual functions 
together to get the overall time-time correlation function.c3 

This allows us to break up the time-time correlation function into the contributions from avalanches of different 
sizes, S. Let G{0\S) be the average time-time correlation function of an avalanche, given that the avalanche is of size 
S. (In contrast, the notation G{9,S) would denote the contribution of all avalanches of size S to G{9). This would 
be weighted by the probability that an avalanche is of size S*, S~'^ .) If we consider equation || at 6* = 0, we see that 
the 9 — Q component of the correlation function is J V(ty'dt, which is proportional to the avalanche energy. Using 
this fact along with the scaling of the energy from equation ^ we find that the time-time correlation function should 
scale as 

G{9\S) = S^-^'''hr..corAOS-^n- (9) 

As shown in figure ^, this scaling works very well over a wide range of avalanche sizes. 

E. The scaling of the power spectrum with S 

The power or energy spectrum describes the amount of energy released in Barkhausen noise at each frequency. It 
can be calculated as the cosine transform of the time-time correlation function. Transforming equation ^ we find that 
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FIG. 4: A collapse of avalanche time-time correlation functions according to equation ||. The data is binned logarithmically to 
get good statistics. 



the scaling of the power spectrum with avalanche size has the form 

/>oo 

E{uj\S) = / cos{uj9)G{9\S)de 



COs{cO0)S^~^'''ftuaccoA0S-''''')d0 

_ q2 r ( l/(7UZ Q\ 

— ^ ^timccorrV'-^ 



(10) 



(We'll call this the power spectrmii, but use the name E[ijj) to remind us to divide by the duration of the experiment 
to change energy to power.) A collapse of the power spectra for different avalanche sizes according to this form is 
shown in figure pi 



F. The power spectrum scaling depends on the large avalanche cutoff 

The power spectrum for all avalanches (the quantity usually measured in experiments) is just the integral 
/ D{S)E{Ljj\S)dS , where D{S) ~ S~'^ is the probability per unit volume of having an avalanche of size S. If we 
substitute equation and D(S) ~ S^'^ into the integral, we generate the naive prediction of 

i?wrong(c^) = f DiS)E{u;\S)dS ^ Cc>-(3-r)/a.. [ y^- f,^,,^{y)dy ^ c^'O-)/-- (11) 



as predicted by Lieneweg et aL.Ill Dahmen and Sethna,Bi@ and Spasojevic et azJll. In the next section [II G we will 
argue that the energy contributed by an avalanche at a fixed frequency to is proportional to S: this immediately 
implies using the scaling form |l^ that /cnorgy(2/) ~ A/y for large y. Thus the integrand in equation |l^ dies off as a 
power law y^^'^ for large y. Hence, it gives the correct answer if r > 2, but for t < 2 the naive indefinite integral in pi] 
diverges and must be replaced by a definite integral. The naive result should work for the integrated avalanche size 
distribution for the short-range modeloCI at the critical point R^- In that case (which we won't study in this paper), 
the corresponding exponent t + al35 — 2.03 in three dimensions (close to two, so we'd expect logarithmic corrections), 
and rising to the mean-field value of 9/4 in six dimensions. For the models we study here, r < 2 and the cutoff at the 
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FIG. 5: A collapse of the power spectra for different avalanche sizes according to equation hol 



largest avalanches in the integral over avalanche sizes changes the form of the pov^fer spectrum. In the section [II H 
after next, we will derive our resulting exponent relation for the power spectrum. 



G. The energy at frequency ui scales linearly with the avalanche size S. 

We shall argue in this section that except at very small frequencies the energy at a fixed frequency lu is proportional 
to S. According to equation |lO|, this implies that the power spectrum at fixed S scales as uj~^^'"^'^ . This means that 
dividing E{lj\S) by Slj^^^'^'^^ should collapse the curves and eliminate the w dependence. The result of this collapse 
for the infinite-range model can be seen in figure ^. 

The fact that the energy at frequency lu scales linearly with S follows from two assumptions. First, each spin in 
the avalanche contributes equally to E{uj). Second, the contribution of a given spin is independent of the size of the 
avalanche it is in. (Of course, a few spins at the beginning of the avalanche might differ, so long as the fraction of such 
spins vanishes for large avalanches.) One way in which this could be true would be if the following two hypotheses 
were true. First, each spin contributes to E[lij) only through its correlations with physically nearby spins. Second, the 
local growth of an avalanche does not reflect the overall size of the avalanche. The first hypothesis seems likely in the 
models described in this paper because the times of physically distant spin flips are likely to be randomly distributed 
in time and contribute incoherently to the power spectrum. Only nearby spin flips will be correlated in time and 
contribute coherently. The second hypothesis is also likely to be true because the avalanches are occurring at a critical 
point. Every part of the avalanche is always on the verge of stopping regardless of the size of the avalanche. 

We can check these hypotheses by breaking up the power spectrum into the contributions from pairs of spins at 
different radii r. The time-time correlation function for avalanches of size S can be written in terms of individual 
pairs of spins as 



G{e\s) 



hi 



(12) 



where ti is the time at which spin i flips. From this form, we can use an additional delta function to pull out the 
contribution due to pairs of spins separated by a distance r: 



G{e,r\S) = (Y,Sit,-U-9)Si\r, 



(13) 
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FIG. 6: A collapse of the power spectra for several avalanche sizes S. The collapse is performed by dividing out the linear S 
dependence. The curves are made flat by dividing out the simple power law a;^'''^''^. The bj axis has not been rescaled, so the 
high uj cutoffs collapse together and the small ci; cutoffs do not. 



where ()i implies an averaging over all values of i. Using this definition of the function G(B^ ''|<5'), we can rewrite the 
time-time correlation function as 

G{B\S) = S J G{e,r\S)dr. (14) 

Now, we can use equation |lj to calculate the contribution of spins separated by a distance r to the power spectrum 
E{uj\S) at a frequency uj. Taking the cosine transform of equation we find that 



E{lu\S) = S J J co8{ioe)G{e,r\S)drde 

= S I E{uj,r\S)dr, (15) 



where the function E{uj, r \S) is defined by this equation. Notice that E{uj, r \S) must have a cutoff at the largest r 
present in an avalanche of size S. We can see from equation |l^ that in order for E(uj\S) to be proportional to S, 
the integral must be independent of this S dependent cutoff. This is a more precise statement of the hypothesis that 
only correlations between nearby spins contribute to the power spectrum. It is also necessary that except for extreme 
values of uj and r, E(uj,r \S) must be independent of S. Combined with the first condition, this corresponds to the 
hypothesis that the local growth of the avalanche should not refiect the overall size of the avalanche. 

Figure ^ shows E{uj,r \S) is decaying approximately as 1/r. Figure || shows that it is independent of S and that 
the decay is oscillating about zero. Neighboring radial shells contribute with opposite sign. Hence the integral |l^ 
appears to be conditionally convergent at large distances: each spin contributes to E{uj) only through its correlations 
with nearby spins, and E{uj\S) ~ S. 



H. Integrating the power spectrum over avalanche sizes 



In section IIIF we found that we need to understand the cutoff at large avalanche sizes. In our simulations there is 
a typical largest avalanche size S* which is primarily due to finite size effects: S* ^ L^l"'^ . In experiments, the typical 
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FIG. 7: This graph shows the function \E{Lj,r \S)\ for a range of values of r at 5 = 32500. The function E{uj,r\S) not only 
decays with uj, but also oscillates about zero. To better compare the amplitudes of the curves at different r, we took the 
absolute value and performed a rumiing average over u, averaging out the oscillations. The horizontal lines show how the 
positions of these curves should scale if they went as 1/r. Notice that while for small r the amplitudes drop slower than 1/r, 
for larger r the amplitudes drop off approximately as 1/r. 



10" 



10" 



10 



-4 



10" 



10" 



S=16250 (r=6) 
S=32500 (r=6) 








0.1 



0.2 



0.3 



0.4 



0.5 



CO 



FIG. 8: This graph shows r \S)\ at r = 6 for two sizes of avalanches. No averaging has been performed, so the oscillations 

are visible (the function changes sign at each dip). For larger r, the oscillations are much faster. Notice that the curves are 
nearly identical for the two sizes — E{u!,r \S) has no significant S dependence at small r. 
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largest avalanche size could be determined by the finite experimental duration or from the demagnetizing forces. (In 
some cases, demagnetizing forces also contribute to the cutoff in our simulations.) These effects would cut off the 
probability of getting large avalanches. 

The contribution to E{lu) from an avalanche of size S scales linearly in S, for large avalanches: each spin contributes 
the same amount. The cutoff at large avalanche sizes dominates the scaling of E{uj) precisely because most spins 
are in the largest avalanches. The distribution of avalanche sizes D(S, S*) cannot be a simple power law cut off at 

S* for T < 2, because we know / SD{S,S*)dS = 1; J" S^-'-dS depends on S* and diver ges as 5* — > oo for t < 2. 
Hence the overall amplitude of D must decreases as S* gets bigger. If we choose a scaling form for the avalanche size 
distribution with an overall amplitude which depends upon S* , we can fix this problem: 

D{S, S*)dS = {S*r^UUS/S*) (16) 

fsi.eiy)-y-^ y^O (17) 

J SD{S, S*)dS = J {S*)-^SU,c{S/S*)dS = J yUUy)dy = 1 (18) 

where the last equation provides a condition on the scaling function /size. 

We are now ready to integrate the power spectrum over S*, finding the form E[uj) corresponding to equation |ll| 
but valid for r < 2. We discovered in the last sections that /encrgy(j/) ^ -^/u- We control the integral by adding and 
subtracting A/y: 

E{uj)= j D{S, S*)E{uj\S)dS 

[{S*)-^UUSIS*)][S^f.n..sy{^^''''''S)]dS 

{S/S*fUUS/S*)foncrsy{^'/"'''S)dS 

{S/S*fUUS/S*)A/{u;^/^'''S)dS 



+ I iS/S*ffsUS/S*)[fener,yiio'/^'''S) - A/ (cj' S)]dS (19) 



In the first integral, we set z = S/S*. In the second integral, we set y = uj^/'^'^^S. Also, the second integral now 
converges, so we can substitute fsizeiS/S*) (S/S*)^'^: 

E{lo) = cc^-i/-- I Az^UU^.)dz + (5*)-2^-(3-r)/a.. f y^-r^j^^^^^^^y) _ Aly]dy (20) 



(t < 2) (21) 

The exponent in the second term of equation (3 — t)/ avz, is -Ewrong of equation given by ignoring the cutoff 
at large avalanche sizes. However, for t < 2 the first term will dominate over the second term both for large system 
sizes and for large oj. Only for r > 2 will the second term dominate. For all of the models we are considering, r < 2. 
(See table |.) In fact, the exponent (3 — T)/avz disagrees badly with the results observed in both simulations and 
experiments, and the exponent l/avz agrees very well. In mean field theory, we can actually derive rigorously that the 
power spectrum exponent is 2 (see appendix |A[). This is in perfect agreement with the exponent l/avz, and complete 
disagreement with the exponent (3 — T)/avz, which would be 3. A plot of the power law for the infinite-range model 
can be seen in figure ^. A comparison of the two possible exponents with simulations for each of the three models can 
be seen in table ||. 



IV. HOW UNIVERSAL IS THE EXPONENT 1/auzl 



So far in this paper, we have only discussed these results for the infinite range model in three dimensions. However, 
we have also checked these results in several other situations where both r and avz take on a range of values. (In 
all cases, r < 2.) We have checked the short-range model in 3 and 4 dimensions, the infinite-range model in 3 and 4 
dimensions, and mean field theory. In all cases, the linear scahng with S seems exact, and both E[ijo\S) and E{uj) scale 
as 0;"^/°''^^ to within simulational precision. In all cases, the exponent (3 — t)I avz is completely inconsistent with the 
observed results. The results of mean field theory and the short-range model at the three dimensional critical point 
can be seen in figures |l^ and |ll]. The results for the four dimensional short-range model and the four dimensional 
infinite range model were also completely consistent with a Xjavz scaling exponent. 
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FIG. 9: The power spectrum for the infinite range model. The dashed line is not a fit. It is a power law with an exponent 
Xjavz. 





Ijavz (3 


~ r) javz 


Simulated 


Short-range 


1.72 


2.41 


1.70 


Infinite-range 


1.72 


3.00 


1.70 


Dipole 


2.00 


3.00 


2.00 



TABLE II: This table compares values of the power spectrum exponent measured in simulations with the values predicted by 
the exponents Xjavz and (3 — T)/auz. Notice that for all three models 1/auz is very close to the simulated exponent, and 
(3 — T^javz is in complete disagreement. 



V. DISCUSSION 

There have been several previous predictions of the power spectrum exponent. In her Ph.D. thesis, Karin DahmenE 
did a calculation similar to the one done in this paper, but she ignored the problems with the integral and came up 
with tiie exponent (3 — t)/ avz, which we have shown should hold only if t > 2. This exporvB«|t-was published by ohs 
groupO without the derivation, and was compared to experimental work by Cote and MeiseOEa and Bertotti et alH 
The correct exponent form derived here makes the agreement between the short-range model and these experiments 
substantially better: they quote an exponent of "around 2", our former (wrong) prediction was 2.46 ± 0.17, and the 
correct prediction is 1.75 ± 0.25. 

Presumahlji the earliest derivation of (3 — T^javz (we thank the referee for pointing this out), is b« Lieneweg and 
Gross-NobisE3 in 1972: their exponent 1 -I- 2e — (5 translates precisely to (3 — T)/avz in our notation. t3 They use this 
to analyze their unannealed data, for which they measure r = 2.1 (hence making their derived exponent appropriate); 
their annealed wires show values of r = 1.73, and indeed their power spectrum shows a much steeper falloff, quite 
close to our predicted value \/uvz which thw, measure to be 1.63. 

Spasojevic, Bukvic, Milosevic and Stanleyli^ also came up with the same (wrong) form based on arguments about 
the average pulse shape. They measured an experimental power spectrum exponent of 1.6. For their experimentally 
determined values of r = 1.77 and avz = 0.662, they found that (3 — T)/avz — 1.86 fit their data well, but the value 
of (3 — T)/avz 2.46 quoted by Dahmen et. aLB was much too large. This lead them them to disregard the plain 
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FIG. 10: The mean field model. Notice that the power spectrum for avalanches of size S collapses when divided by S, and 
that the power spectrum for avalanches of size S and the overall power spectrum both have the same power law as the average 
avalanche size as a function of the inverse avalanche duration. Sizes 96, 1536, 24576, and 196608 are shown in the collapse of 
P{uj\S)/ S. The average avalanche size as a function of the inverse avalanche duration is the definition of the exponent 1/auz. 
Power laws of l/ai/z and (3 — T)/auz are also shown for comparison. Notice that (3 — T)/auz is completely incompatible with 
the results. 




FIG. 11: The short-range model in 3 dimensions. The plots and the sizes are the same as in figure y_0| 
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old critical model with domain nucleation as a possible explanation of Barkhausen noise in their system. However, 
note that their experimental value of X/avz — 1.51 is as close to their measured power spectrum exponent as their 
value of (3 — t)/(7vz ~ 1.86. The value of 1/avz — 1.75 ± 0.25 quoted by Dahmen et al. is even closer to the 
experimentally observed power spectrum exponent. (The value of t = 1.6 predicted by Dahmen et al. is also closer to 
the experimentally measured r — 1.77 than any of the other models: front propagation has an exponent r = 1.3, and 
mean-field theory has an exponent r = 1.5.) One should note, however, that the average avalanche shapes measured 
by Spasojevic et al. disagree with those prec^ted by the short-range model (figure 

In 1989, Jensen, Christensen and FogedbyEZI published a different calculation of the power spectrum exponent for 
sandpile models, which has been cited many times as an explanation for 1/w^ scaling of the power spectrum in 
sand-piles, Barkhausen noise, and other systems. They made two major assumptions in their derivation. First, they 
assumed that the avalanche shape could be approximated by a box function: V{t) ~ S/T for all t < T. For our models, 
this assumption turns out to be valid for calculating the average avalanche energy, but is not valid for determining the 
overall scaling of the time-time correlation function. Second, they assumed that one of their scaling functions, which 
was related to the time-time correlation function, had the simple scal ing form G{T) ~ T" exp(— T/Tb). As can be seen 
from the exact mean field time-time correlation function in equation A6, this is not at all a safe assumption. In terms 
of our exponents t and ai^z, their prediction was that the power spectrum would have the form E{uj) ^ uji^^'^)/'^'^^ 
for r > 3 — 2ai'z and E{uj) ^ uj^'^ for r < 3 — 2(ji'z. This is the same as our result in the case where avz — 1/2. 
For the particular model they were studying, their assumptions seem to be correct, but in the more general case, the 
exponent (3 — T)/avz should become invalid for r < 2, rather than r < 3 — avz. Also, for t < 2, exponents other 
than 2 are possible, depending on the particular dynamics. For the class of models described in this paper, the correct 
exponent for. r < 2 is l/avz. 

NarayanES derives a power-spectrum exponent for fcont propagation, considering the dynamic height- height correla- 
tion function for the growing interface. His exponentiHa, {2(^ + d— l)/z — 1 involves the static height-height correlation 
exponent C: {h{x,t)h' {x' ,t)) ~ jx — x'p'', along with the dimension d. This^rediction is rather different from the two 
discussed in this paper. Narayan and Fisher's 5 — d = e-expansion resultsB2l show ( = 2/3 and z = 2 — 2e/9 + O(e^) 
yielding a predicted power-spectrum exponent of 8/7 -I- O(e^) in three dimensions, quite different from the value of 
around 1.7 that we measure. Narayan's derivation does not separate the high-frequency power law 'Pav(t^) which 
we study from the low-frequency power law 7'corr('^) reflecting the correlations between avalanches (discussed in the 
introduction): his analysis is probably describing the asymptotic behavior of the latter. 

There are other specific models for which the power spectrum exponent has been calculated. For example, Bak, 
Tang and WiesenfeldH calculated the power spectrum exponent for their sandpile models, and—Diirin, Bertotti, and 
Magncj calculate the form and asymptotic power law for the power-spectrum of the ABBMMEj model. Without 
further investigation, we can't expect our arguments for an exponent of l/avz to hold for these and other models. 
However we do expect that in any avalanche based model, the power spectrum exponent for r > 2 will be (3 — r)/ avz, 
and another exponent will dominate for r < 2. Whenever the arguments in section III G hold, for r < 2 the exponent 
l/avz, relating the avalanche size to avalanche duration, should dominate. 
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APPENDIX A: THE MEAN FIELD POWER SPECTRUM 



In mean field theory, we can calculate the power spectrum exactly. The Hamiltonian in mean field theory is the 
Hamiltonian in equation ^ without nearest neighbor or dipole terms: 

H = -Y2iH + J + h^)s^. (Al) 

i 

When a spin flips with an external field of H, all spins with random fields between —{H+JM) and —{H+J{M+2/N)) 
will flip. Therefore, each spin has a probability of ^p{—H — JM) of flipping, where p{h) is the probability distribution 
of the random fields. On average, 2Jp{—H — JM) spins will be flipped. If 2J p{—H — JM) > 1, then the avalanche 
will tend to grow indefinitely, and there will be an infinite avalanche. If 2J p{—H — JM) < 1, then the avalanche 
will quickly die out, and all avalanches will be small. If 2Jp{—H — JM) — 1, then avalanches will always be finely 



14 



balanced between continuing and dying, and there will be a critical distribution of avalanches. If the random field 
distribution p{h) has a maximum value of 1/2 J, then there will be a critical distribution of avalanches at the value of 
H = He where p{h) is a maximum. 

Let us calculate the power spectrum for the critical system with p{h) ~ 1/2 J. To begin with, we will calculate the 
probability distribution of time-series rti,n2, . . . ,noo- We know that in shell zero, exactly hq — I spins will flip. At 
the critical point, where each spin on average causes one more spin to flip, shell one will have a Poisson distribution 

with mean one: P{ni) = Shell i will have a Poisson distribution with mean ni_i: P{ni) — Therefore, 

the probability distribution for the entire time series will be 



2 °° e"*"in"' 

P(l, n„n„...,n^)^—J[ (A2) 



Now, from equation A2, we can calculate the average time-time correlation function 

oo 

G{e) = (A3) 

i=0 

where 

{n^Ui+g) ^ ^ nin^+gP{ni,...,noo)- (A4) 

{nl,...,n^}=0 

To simplify equation |A^ , we need to use several properties of the Poisson distribution: J2^=a ^ „f = 1, 

X^^o*^^ — r' ~ ^n=o T^r" ~ X + x^. Using the first rule repeatedly, we can simplify equation A_4 

to 

OO 

{niUi+g) ^ ^ nifli+gPim, Tli+g). 

{ni,...,ni+8}=0 

Then, using the second rule repeatedly, we can further simplify to 

oo 

{niUi+g) = ^ nfP{ni,...,ni). 

{ni ,...,ni}— 

Now, applying the second and third rules repeatedly, we can simplify to a single sum: 

{i — 1)77,1 + ni 



(jliTli+g) = ^ 



eni ! 

rii=0 ^ 



= i + l. (A5) 

Notice that the correlation between two times is proportional only to the first time, and not the separation between 
the times. 



Now, we can find the value of the time-time correlation function G{9) by summing the result of equation A5 



Because the time-time correlation function as defined is proportional to the square of the total time, we must cut off 



the summation at some maximum time T to get a finite result. Summing according to equation A3, we find 

i=0 

_ {T -e){T-e + i) 
~ 2 

2222 ^ ' 

This shape for a cutoff time of T = 1000 is shown in figure |l^. Notice that the exact shape in mean-field theory is 
very similar to the experimentally measured correlation functions in three dimensions. 
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FIG. 12: The exact time-time correlation for mean field theory with a maximum time cutoff of T = 1000. 



To get the exact form of the mean-field power spectrum, we just take the cosine transform of equation A6 

E{oj) = / cos{ujd)G{0)de 
Jo 



Tlj- 



(1 - cos(rw)) - uj-^ sin(rw). 



(A7) 



The dominant term in this equation is Tlu^'^. This is the same exponent as predicted by the general scaling arguments: 
-^j^ = 2. The general scaling arguments also predicted that a term smaller by a factor of T and with the exponent 
(3 — t)I<jvz = 3 would be subtracted off, and it is, but multiplied by a factor of sm^TtS). However, the sin(rti;) turns 
out to simplify things even more. Because the correlation function was actually discrete, we should consider a discrete 
power spectrum, where the frequencies are multiples of = This means that cos(Ta;) = 1 and sin(Tci;) = for 
all u) in the discrete spectrum. Because of this, all terms except the Tuj~'^ term drop out. If we divide by T to get 
the power spectrum, we find 



(A8) 
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